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The equation of the Van der Pol oscillator, being characterized by a dissipative term, is non- 
Lagrangian. Appending an additional degree of freedom we bring the equation in the frame of 
action principle and thus introduce a one-way coupled system. As with the Van der Pol oscillator, 
the coupled system also involves only one parameter that controls the dynamics. The response 
system is described by a linear differential equation coupled nonlinearly to the drive system. In the 
linear approximation the equations of our coupled system coincide with those of the Bateman dual 
system (a pair of damped and anti-damped harmonic oscillators). The critical point of damped and 
anti-damped oscillators are stable and unstable for all physical values of the frictional coefficient fi. 
Contrarily, the critical points of the drive- (Van der Pol) and response systems depend crucially on 
the values of /i. These points are unstable for fi > while the critical point of the drive system is 
stable and that of the response system is unstable for /i < 0. The one-way coupled system exhibits 
bifurcations which are different from those of the uncoupled Van der Pol oscillator. Our system is 
chaotic and we observe phase synchronization in the regime of dynamic chaos only for small values 
of fl. 

PACS numbers: 05.45.-a, 05.45.Xt 
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I. INTRODUCTION 

A great variety of periodic processes, ranging from re- 
occurrence of economic crisis to beating of the human 
heart, takes place in such a way that external energy 
is supplied to the system over a part of the period and 
dissipated within the system in another part of the pe- 
riod. Periodic processes of this type have been qualified 
as relaxation oscillations. In fact, the relaxation oscil- 
lator constitutes a self-sustained system with a built in 
mechanism to pump energy to oscillations that are too 
small and damp those oscillations that grow too large. It 
appears that Van der Pol [l| was the first to model relax- 
ation oscillations in a laboratory experiment by consid- 
ering an electric circuit which consists of a battery con- 
nected in series with a resistor and a capacitor. A diode 
valve which acts as a nonlinear device was also connected 
across the capacitor. 

Because of their applicative relevance, the driven and 
coupled Van der Pol oscillators have been extensively 
studied in the literature 0. These oscillators display 
physical and mathematical properties which are radically 
different from those of an isolated Van der Pol oscilla- 
tor. In the present work we are interested to construct a 
model for the one-way coupled system Q in which Van 
der Pol oscillator is a drive. Initiated by Chua [4|], studies 
in the dynamical properties of one-way coupled systems 
have become a subject of considerable interest because 
such systems exhibit bifurcation, chaos and synchroniza- 
tion. Thus we can expect that studies in our newly pro- 
posed system may also lead to similar physically signifi- 
cant results. 

In section 2 we note that Van der Pol oscillator is not 
invariant under time reversal and thus represents a physi- 



cally incomplete system Q in the sense of the Hamilton's 
variational principle. We follow Bateman [|| and double 
the degree of freedom to bring in an associated equation 
such that the system of equations follows from the action 
principle. Curiously enough, this automatically provides 
us with a one-way coupled system in which Van der Pol 
oscillator is a drive. In section 3 we write the equations 
for coupled Van der Pol system as a system of first or- 
der differential equations that provide a basis to study 
the associated dynamics. For an uncoupled Van der Pol 
oscillator the stability of critical points is well known. 
It has a limit cycle attractor and bifurcation structure 
dependent on the parameter of the oscillator. We shall 
study the stability of our one-way coupled system with 
particular emphasis on the nature of the attractor and 
bifurcation that lead to chaos. One of the remarkable 
features of coupled nonlinear oscillators is synchroniza- 
tion which characterizes the emergence of coherent mo- 
tion among the constituent oscillators. Keeping this in 
view we shall also analyze the stability and optimization 
of the synchronization process which is likely to occur 
in our system. Finally, in section 4 we shall make some 
concluding remarks. 



II. ONE-WAY COUPLED SYSTEM : ACTION 
PRINCIPLE AND LINEAR ANALOG 



.1] 



The equation of the Van der Pol oscillator is given by 



G(x, x, x) = x — — x 2 )x + x = 0, 



(1) 



where overdots denote differentiation with respect to 
time and /z is a parameter of the system. Clearly, this 
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oscillator is driven by a velocity dependent force and pro- 
vides a typical example of the dissipative system. More- 
over, the equation violates the time reversal symmetry. 
Such equations do not follow from the action principle 
to have a Lagrangian representation. It is well known 
that a dynamical equation can be derived by using the 
Hamilton's variational principle if it satisfies the criterion 



dG(x,x,x) d dG(x,x,x) . . 

dx dt dx 

It is easy to see that (1) invalidates the condition in 
(2) such that the equation of the Van der Pol oscilla- 
tor is non-Lagrangian. A non-Lagrangian system can be 
brought within the framework of variational principle by 
doubling the degrees of freedom of the equation of motion 
Q. To illustrate this consider the time reversal symme- 
try breaking equation of the linearly damped harmonic 
oscillator 



similar to that in (6). As opposed to the uncoupled equa- 
tions (3) and (6) , we have now a coupled system in which 
the Van der Pol oscillator drives the associated response 
system through nonlinear coupling. All contemporary 
physical theories use Hamilton's variational principle to 
decode the laws of nature and the one-way coupled sys- 
tem obtained by us follows from the constraint imposed 
by this principle. Looking from this point of view the 
equations in (1) and (8) which represent a one-way cou- 
pled system appear to be a natural selection. This is not 
true for other widely cited examples of one-way coupled 
systems. In this context we note that the response sys- 
tem in (8) is governed by a linear differential equation 
coupled nonlinear ly to the Van der Pol oscillator. In- 
terestingly, if this coupling were absent and if we would 
linearize the equation of the Van der Pol osillator, we 
will get back the equation for the Bateman dual system. 
Thus the Bateman system can be considered as a linear 
analog of our one-way coupled system. 



v + m + v = o. 



(3) 



The quantity fi now stands for the dissipation constant 
of the medium in which the oscillator is embedded. Let 
us write a Lagrangian 



L = x(y + ny + y), 



(4) 



characterized by two degrees of freedom, namely x and y. 
Clearly, the Euler-Lagrange equation in x gives (3). To 
obtain the equation of motion for the x variable consider 
the Euler-Lagrange equation [8j 



f_ fdL 

dt 2 \ dy 



d fdL 
Jt \~dy~ 



dL 
dy 



0. 



in the y variable. From (4) and (5) we get 
x — fix + x = 0. 



(5) 



(6) 



Equation (6) models a physical system which absorbs 
energy dissipated in the damped harmonic oscillator (3). 
The system represented by (3) and (6) is often called 
dual system of Bateman who initiated the idea of bring- 
ing non-Lagrangian systems within the frame of action 
principle. However, it is not immediately clear how the 
transfer of energy takes place between the components of 
the system. In the following we show that the approach 
of Bateman when applied to the equation of a Van der 
Pol oscillator we get a physically transparent one-way 
coupled system with the Van der Pol oscillator as the 
drive system. 

To bring the equation of motion (1) in the frame of the 
action principle let us define a Lagrangian 



L = y(x — — x 2 )x + x), 
so as to write an associated equation 

jj + fJ,(l -x 2 )y + y = 0, 



(7) 



(8) 



III. STABILITY ANALYSIS 

For stability analysis and even for purposes of numeri- 
cal simulations, it is convenient to write (1) and (8) as a 
system of first order differential equations 



xi = x 2 , 



x 2 = ^(1 - x\)x 2 - X\, 



Vi = 2/2, 



and 



y-2 



- x i)V2 - Vi- 



(9a) 



(9b) 



(10a) 



(10b) 



In writing (9) and (10), we have used x\ for x and y\ for 
y. Clearly, x\ — x 2 — yi — yi — represents the fixed 
point of (9) and (10). Linearizing the system about the 
fixed point , say, £ and rj we have 



Cl - C2, 



C2 = K2 - Ci, 



and 



<12 



(11a) 



(lib) 



(11c) 



(lid) 



In writing (11), we have used the same notational conven- 
tion as that adopted for x and y. The Jacobian matrix 
for (11) given by 
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J = 



has the eigenvalues 








1 
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1 










-1 





(12) 




Ai = 1+2^-4, (12a) 

M = \-\^^l, (12b) 

A 3 = -^ + ^V / ^^, (12c) 

and 

^ = ~\~\^^- (12d) 

The nature of the eigenvalues gives a broad classification 
scheme for the equilibrium points regarding their stabily 
and instability. We shall use the method for such classi- 
fication wherever necessary. 

If the nonlinear terms in (1) and (8) were absent the 
drive and response systems will go over to the Bateman 
equations (6) and (3). Thus it will be interesting to com- 
pare the phase diagrams of (1) and (6) and of (8) and (3) 
in order to study the points of contrast or of similarity for 
the phase-space evolution of the one-way coupled system 
and Bateman dual system. 

To draw the phase diagrams of Bateman system we 
rewrite (3) and (6) as a system of first order differential 
equation 

2/1=2/2, (13a) 
2/2 = -M2/2 - 2/i, (13b) 

and 

±i = x 2 , (14a) 

±2 = [iX2 — x\. (14b) 

As with (9) and (10), we used yi — y and x\ = x to 
write the equations in (14) and (15). We solved this 
set of equations using the initial conditions yi(0) = 1, 
y 2 (0) = and .ti(0) = 0.2, x 2 (0) = 0. We have chosen to 
work with fi = 0.2. The phase diagrams of the damped 
harmonic oscillator (3) is given in Fig. 1(a). Here we see 
that the phase trajectory starting from the point (1,0) 
spirals in and tends to meet the origin such that the fixed 
point is a stable focus.We have presented in Fig. 1(6) a 
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FIG. 1: Phase diagrams of Bateman dual system. (o) Damped 
harmonic oscillator and (b) The associated system. 



similar phase diagram for the associated equation in (6). 
The phase trajectory starts from (0.2, 0) and spirals out. 
Thus the fixed point is an unstable focus. We shall now 
make similar comparison between the phase diagrams of 
(1) and (8) .We solved the system of equations (9) and 
(10) with the same initial conditions as used for (13) and 
(14). Fig. 2(a) gives the phase diagram for the Van der 
Pol oscillator. As expected the phase trajectory shows 
a stable limit cycle with origin as an unstable focus for 
/i = 0.2. The phase diagram (Fig. 2(6)) of the response 
system, which in conjunction with the Van der Pol os- 
cillator forms a one-way coupled system, exihibits some 
interesting features. For example, starting from (1, 0) the 
phase trajectory tends to form a small amplitude limit 
cycle. Then all of a sudden it spirals out. 

Since in the absence of the nonlinear term, (1) concides 
with (6), we will compare the phase trajectory in Fig. 2(a) 
with that in Fig. 1(6). Looking closely into these figures 
we see that the phase trajectory in Fig. 1(6) spirals out 
without forming a limit cycle while the phase trajectory 
in Fig. 2(a) has a stable limit cycle. Thus the forma- 
tion of limit cycle is a typical feature of sustained oscilla- 
tion. Again comparison between the phase trajectories of 
Fig. 2 (6) and Fig. 1(a) leads us to more surprising points of 
contrast because the phase trajectory in Fig. 1(a) simply 
spirals in rather than showing any complicated behavior 
as found for the phase trajectory in Fig. 2(6). Ideally, one 
would expect some similarities between the phase tra- 
jectories of Fig. 2(6) and Fig. 1(a). This is, however, not 
true. The phase-space evolution of the drive system is 
radically different from that of the damped harmonic os- 
cillator because the energy transfer between the systems 
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FIG. 2: Phase diagrams of the one-way coupled system, (a) 
Drive system (Van der Pol oscillator) and (b) Response system 
(Equation (8)) 



represented by (1) and (8) is explicit while the equations 
of the Bateman dual system are uncoupled. 
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FIG. 3: Maximum Lyapunov exponent versus parameter of 
the system 



We have drawn the phase diagrams of our one-way 
coupled system for fi = 0.2 which satisfies the inequal- 
ity < ji < 2. If /i is constrained by this inequality 
the uncoupled Van der Pol oscillator will not exhibit any 
bifurcation. However, there are three critical values of 
ji for which certain qualitative changes, namely, bifur- 
cations occur in the dynamics of (1). For example at 
// = —2 there occurs a bifurcation in which the criti- 
cal point changes from a stable node to a stable focus. 
The radicals in Ai and A2 (eigenvalues of the equation 



for the uncoupled Van der Pol oscillator) do not depend 
on the sign of \i. When the real parts of the eigenvalues 
go through zero {p, = 0) with non-zero complex parts we 
have Hopf bifurcation. At fi = 2 a bifurcation from an 
unstable focus to an unstable node occurs. 

Consider a similar situation for the coupled system. 
The response system represented by (8) has also three 
bifurcation points. For fi > — 2 the phase trajectory di- 
verges from the critical point such that there is no attrac- 
tor in the neighborhood of the equilibrium point. As we 
increase the value of \i and choose to work in the region 
— 2 < ji < the equilibrium point becomes an unstable 
focus similar to that in Fig. 1(6). Thus at yU = —2, a bi- 
furcation takes place in which a diverged state changes 
to an unstable focus. If we still increase the value of /i 
through zero (0 < /i < 2) there appears a fold limit cycle 
(Fig. 2 (6)). This small amplitude limit cycle is neither 
stable nor unstable. As noted earlier the phase trajec- 
tory first spirals in, tends to form a limit cycle, then all 
of a sudden spirals out. This implies that at ji = a 
bifurcation occurs from a unstable focus to a fold limit 
cycle. For \i > 2 the fold limit cycle disappears and the 
phase trajectory again diverges from its critical point. So 
at n = 2, a bifurcation occurs in which a fold limit cycle 
vanishes and the state of the system diverges. 

In the above we have seen that both uncoupled Van 
der Pol oscillator and response system driven by it ex- 
hibit bifurcations at /i = —2, and 2 although the type 
of bifurcation is different in each case. Since the appear- 
ance of bifurcation in a dynamical system signals onset 
of chaos, it will be interesting to examine which of these 
system is chaotic for — 2 < [i < 2. 
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FIG. 4: x(t) (thin line) and y(t) (thick line) as a function of 
t 



The non-wandering set is a set of points in the phase 
space for which all orbits starting from a point of the 
set come arbitrarily close and arbitrarily often to any 
point in the set. The different types of non-wandering 
sets include fixed point, limit cycle, quasi-periodic orbits 
and chaotic orbits. A non- wandering set can be either 
stable/unstable as one changes the parameter of the sys- 
tem. In particular, such a set is said to be Lyapunov 
stable if every orbit starting in its neighborhood remains 
in its neighborhood else it is unstable. The linear sta- 
bility analysis of the system leads, in a rather natural 
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way, to the concept of Lyapunov stability. Let X(t) be 
an orbit satisfying 



X(t) = f(X,t). 



(15) 



A solution X(i) of (15) will be asymptotically stable if 
any infinitesimally small perturbation 8X(t) decays. In 
order to determine the stability of any solution Xo(t) of 
(15) assume 



X(t) = X (t) + 6X(t), 



(16) 



as one of the purturbed solution of (15). Using (16) in 
(15) we expand the function f(X,t) in Taylor series to 
get 



d_ 



AX 



dX' 



AX. 



(17) 



The linearized equation of motion for AX in (17) is jus- 
tified as long as the orbit is in the neighborhood of X$. 
At the fixed point = a constant (say A). Thus 



AX 



AX. 



(18) 



where, A s and X s are eigenvalue and eigenvector of the 
Jacobian matrix of the linearized system and s is the di- 
mension of the phase space. There are as many Lyapunov 
exponents as dimensions of the phase space. 

The Lyapunov exponent A s can be calculated for each 
dimension. When talking about a single exponent one 
is normally referring to the largest. If the Lyapunov ex- 
ponent is positive the system is chaotic and unstable. 
In this case the nearby points will diverge irrespective of 
how close they are. Thus the magnitude of the Lyapunov 
exponent is a measure of chaosness in the system. 



TABLE I: Signs of X„ 
A 4 



in different ranges of the parameter 



SI. No. 


Range of li 


Sign Of Xmax 


Nature of the system 


1. 


li < -1.5 


positive 


chaotic 


2. 


li = -1.5 


zero 


conservative 


3. 


-1.5 < a* < 


negative 


dissipative 


4. 


li = 


zero 


conservative 


5. 


< (j. < 0.92 


positive 


chaotic 


6. 


H = 0.92 


zero 


conservative 


7. 


0.92 < n < 1.25 


negative 


dissipative 


8. 


li = 1.25 


zero 


conservative 


9. 


H > 1.25 


positive 


chaotic 



In Fig. 3 we plot the Lyapunov exponent, \rnax , versus 
[i, the only parameter of our system. The curve for \ max 
clearly shows that the exponent takes up positive, zero 
and negative values such that the corresponding motion 



of the system is chaotic, neutrally stable (conservative) 
and dissipative. Table 1 gives the signs of \ max in differ- 
ent ranges of the parameter. This helps us identify the 
nature of the system as a function of li. Fig. 2(6) gives the 
phase diagram of the response system for \i = 0.2. This 
li value satisfies the inequality < [i < 0.92(S1. No. 5). 
Thus the curve in this figure gives the phase space evo- 
lution of the system when it is chaotic. 

The response system is synchronized with the drive 
counterpart when the Lyapunov exponent changes sign. 
From Fig. 3 we observe such changes for li values rang- 
ing from —1.5 to 1.25. In Fig. 4 we plot x(t) (thin line) 
and y(t) (thick line), the solutions of the equations for 
the drive and response systems, as a function of time t. 
This Figure shows that during time evolution y(t) is only 
partially synchronized with x(t) in the sense that the am- 
plitudes of y(t) and x(t) remain unsynchronized but their 
phases appear to evolve in synchrony. Synchronization of 
this type, often called the phase synchronization, occurs 
when the coupled oscillators are not identical |9j. From 
(1) and (8), we see that the coupled oscillators of our one- 
way system are nonidentical. As a result the drive and 
response systems are found to be partially synchronized. 



IV. CONCLUSION 

In point mechanics and classical field theory there are 
many physical systems for which the equations of motion 
do not follow the action principle. One of the nonmax- 
imal way to bring such equations in the framework of 
variational principle consists in doubling the degrees of 
freedom of the system. We followed this route to make 
the Van der Pol equation variationally self-adjoint and 
thereby construct a one-way coupled system. The sys- 
tem constructed by us involves two non-identical oscilla- 
tors. Traditionally, the drive-response scenario of chaos 
is studied by using one-way coupled systems composed of 
identical oscillators. The system introduced by us may 
play a role in biology since one-way entrainment occurs 
in heart. A pacemaker entrains the nerve cells of heart 
to prevent heart attacks. Such an applicative relevance 
appears to be further justified by the historically signif- 
icant work of Van der Pol and Van der Mark who used 
coupled oscillators discovered by them to make the first 
model for an artificial pacemaker. 
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